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Abstract 

Orbital Kondo effect is treated in a model, where additional to the conduc- 
tion band there are localized orbitals close to the Fermi energy. If the hopping 
between the conduction band and the localized heavy orbitals depends on the 
occupation of the atomic orbitals in the conduction band then orbital Kondo 
correlation occurs. The noncommutative nature of the coupling required for 



the Kondo effect is formally due to the form factors associated with the as- 
sisted hopping which in the momentum representation depends on the mo- 
menta of the conduction electrons involved. The leading logarithmic vertex 
corrections are due to the local Coulomb interaction between the electrons 
on the heavy orbital and in the conduction band. The renormalized vertex 
functions are obtained as a solution of a closed set of differential equations 
and they show power behavior. The amplitude of large renormalization is de- 
termined by an infrared cutoff due to finite energy and dispersion of the heavy 
particles. The enhanced assisted hopping rate results in mass enhancement 
and attractive interaction in the conduction band. The superconductivity 
transition temperature calculated is largest for intermediate mass enhance- 
ment, m* /m ~ 2 — 3. For larger mass enhancement the small one particle 
weight {Z) in the Green's function reduces the transition temperature which 
may be characteristic for other models as well. The theory is developed for 
different one-dimensional and square lattice models, but the applicability is 
not limited to them. In the one-dimensional case charge- and spin-density 
susceptibilities are also discussed. Good candidates for the heavy orbital are 
/-bands in the heavy fermionic systems and non-bonding oxygen orbitals in 
high temperature superconductors and different flat bands in the quasi-one 
dimensional organic conductors. 
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I. INTRODUCTION 

Recently the different phase transitions in strongly correlated electronic systems has 
became the most intensively discussed issue in the theory of solids. 

There is a common feature in the theories of these phase transitions, namely logarith- 
mic corrections characteristic for infrared divergences occur in the perturbative expansion. 
That genuine character exists even for the phase transitions in electronic systems with the 
logarithmic expression ln{x/D), where x = max(co',T) is the largest of the energy variables 
uj and the temperature T, furthermore, the high energy cutoff D is the bandwidth of the 
conduction electrons. In particular cases an additional infrared, low energy cutoff appears 
which is an inherent feature of some models. That cutoff smears out the infrared divergences 
mentioned earlier. 

Such logarithmic corrections occur in the theory of a one-dimensional electron gasa and 
also e.g. in two-dimensiona. In the first case that is due to the simple Fermi surface, 
while in the second case either the Fermi surface is of one-dimensional character (nested 
Fermi surface) or that is originated by the corners of the Fermi surface for a nearly half-filled 
electron bandB. These special properties are required to get logarithmic terms in the electron 
hole channel (zero sound channel), on the other hand, it is well known from the theory 
of the superconductivity, that the electron-electron channel with total zero momentum 
(Cooper channel) is always divergent. The appearance of the logarithmic contributions are 
generally controlled by the conservation of momentum. The importance of that conservation 
is essentially reduced in those two-band models of arbitrary dimension, where one of the 
bands is almost dispersionless thus the heavy particle can absorb an arbitrary momentum 
with almost the same energy transfer. If the heavy band is close to the Fermi energy then 
the diagrams containing only one light electron (hole) in the intermediate state contribute by 
logarithmic terms arising from the integration with respect to the energy the light electron 
(hole). The finite, but small energy and dispersion of the heavy band introduce an inherent 
smearing of the logarithmic terms, thus it represents a small infrared cutoff mentioned 



earlier. 

Furthermore, in order to get relevant leading order logarithmic corrections in the absence 
of special features of the Fermi surface, it is not enough to have an almost dispersionless 
branch in the excitation spectrum like localized spin excitations, two-level systemsEl a heavy 
band etc., as for the most simple models the leading two logarithmic vertex corrections 
due to intermediate states with a single light electron and hole accompanying the heavy 
electron cancel each other. Such cancellation occur in the theory of the X-ray absorption 
by a deep electron level, which phenomena was worked out by MahanQ, Nozieres and de 
Dominicistm. The logarithmic character of single loop approximation, however, sustains 
in those cases where the couplings between the light and heavy electrons show non-trivial 
structure exhibiting noncommutative behavior which may be due to the 

(i) spin dependenceQ (Kondo effect) or dependence on the atomic orbital indices of a 
single atomic siteOa, 

(ii) dependence on the momenta of the light electrons, which might appear as a structure 
factor involving two sites in the couplingpll2ryj. 

The present paper is devoted to the second case, where the structure factor is associated 
with the electron assisted hopping (correlated hopping) in a localized orbital pictureE^ry. 
In those models the electron hops between two different sites, but the hopping rate depends 
on the occupation of an other site by electrons. The hopping considered may be between 
two localized heavy electron stateal^ or between a heavy and light electrons mixing the two 
bandsEJtJ. In solids the mixing cases have usually larger amplitudes as the hopping between 
two heavy orbitals of smaller radius are essentially weaker. 

The general form of the electron assisted interaction in the real space contains the fol- 
lowing combination of creation and annihilation operators 

where ni, n2 are sites involved by the hopping which is affected by the occupation of site n 
and a, f3, 7 are band indices. In the following only such models are considered where either 



rii = n OT n2 = n. 

That interaction in the Bloch wave form is associated with a from factors just hke in the 

tight binding approximation. 

There are two simple sources for such assisted hopping processes: 

(i) the hopping between two sites must depend on the occupation of these sites by other 

electrons as additional electron or hole on the site modifies the size of the orbital in the 

real space due to the Coulomb interaction, thus the hopping matrix element is changed as 



we' "^^' 



(ii) in the site representation of the two particle Coulomb interaction which contains 
two creation and two annihilation operators there are always such terms, in which one site 
appears three times and another site only onceEj'tB. These terms represent the off-diagonal 
Coulomb interactiont3. 

The form factor appearing in the Fourier transform of the operator product given by the 
expression (|1.1| ) for e.g. ni = n is 

^'^''^''-''-''^'''^''''''<k.^,M<i^.^.M (1-2) 

where the locations of sites ni and n2 are -Ri and R2, furthermore SR = R2 — Ri- The 
first factor in this expression drops out as it ensures the momentum conservation which 
is of limited relevance where one of the bands is flat. The important form factor which is 
responsible for the noncommutative nature of the interaction is the factor e*'^^^ Considering 
the light electron assisted mixing term between the heavy and light bands one of a and /3 
belongs to the heavy band and all the others to the light band. 

The large renormalizations of the different quantities as the strength of the interaction, 
the mass of the light electrons and the electron-electron interaction relevant for supercon- 
ductivity are due to the Coulomb interaction between the electrons in the light and heavy 
band which screens the charge of the heavy electrons by the light ones. The structure of 
these interaction is not crucial, thus it can be taken as a local on-site interaction with 
strength U, except for those cases where the model itself requires more complicated struc- 



tures. Such interactions itself without assisted hopping does not lead to leading logarithmic 
correct ionsllZltJ. The theories belonging to that general class exhibit the following common 
features: 

(i) The vertices are renormalized and they are power functions of the largest of the 
variables 00/ D, T / D and Eq/D with a low energy infrared cutoff Eq discussed earlier. In 
the models where there is only a structureless screening interaction U between the heavy 
and light electrons the leading term of the exponent is at least quadratic in the interaction. 
In the case of electron assisted band mixing treated in the present work the exponent is 
linear in U as that mixing interaction is coupled to three light electronst3, in contrary to 
the renormalization of the screening interactionMllHl. The coefficient of U is different from 
zero even if the structure factor is absent in the interaction, but its presence can enhance it 
by a factor of two. Such enhancement can be very important. 

(ii) The light electron mass corrections exist always due to the assisted mixing, but the 
large vertex corrections can lead to a large enhancement which may reach several orders of 
magnitude (heavy fermionic behavior )li^. 

(iii) The electron-electron interaction induced by electron assisted mixing, where the 
heavy electron occurs in the intermediate state, is different from zero in the presence of the 
form factorlilHiHHU. The vertex corrections are important and the interaction may depend 
strongly on the momenta of the electrons. Without form factors such interactions are not 
generatedllZrEj. 

(iv) Both the induced electron-electron interaction and the mass enhancement increase 
with the strength of coupling in the weak coupling limit. For stronger coupling the density 
of the light electron increases, but the single particle weights Z in their spectral functions are 
drastically reduced. In the induced electron-electron interaction the square of that weights, 
Z^ occurs, which dominates over the increase in the electron density which is proportional 
to Z^-^. Thus in these models the strength of the superconductivity is the strongest for 
moderate mass enhancement 2-3 and decreasing very fast for larger mass enhancement 
{Z ^ 1). The largest available dimensionless attractive coupling is close to unity. This 
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feature must be quite general for models where the electron-electron interaction is induced 
and competes with the mass enhancements. 

In all of the behaviors listed above the energy and the dispersion of the heavy band play 
determining role by limiting the enhancement of the vertex (infrared cutoff). 

It is interesting to note, that in most of the casescl of logarithmic problems new couplings 
are generated by solving the parquet or renormalization group equations for the vertex 
functions. That is the case for the two dimensional electron gas as welln. In the present 
case the finite number of the generated form factors reduces the problem to a closed set of 
differential equations, and their solution will be presented in analytical form. 

The models where the present theory can be relevant are those, where there is a heavy flat 
band. Those bands can be formed by the atomic orbitals of very small size like 4/ electrons in 
heavy fermionic systems and the mixing is between s— and /—electrons. Another possibility 
is the case of non-bonding orbitals in complicated crystal structuresEJ, where these orbitals 
are only weakly hybridized with the conducting electrons, e.g. non-bonding tt— orbital of 
oxygen in the high temperature superconductors or some appropriate orbitals in the quasi 
one-dimensional conducting molecular crystals where there might be several flat bands. The 
present paper is not devoted to study particular special cases, but to provide several methods 
which can be applied to different concrete cases. 

The paper is organized as follows: In chapter II. we are presenting the Hamiltonian of a 
general model with features described above. In chapter III. we are calculating the vertex 
corrections for the general model. In chapters IV. and V. the calculations of the self-energy 
and the electron-electron interaction are presented and in the next chapter (VI.) the su- 
perconducting transition temperature is determined. That discussion of the relation of the 
mass enhancement to the generated electron-electron interaction contains quite general con- 
siderations which are valid much beyond the model treated here. The formalism developed 
in the previous chapters is applied to special cases: to the one dimensional electron gas 
(chapter VII., see also Ref. [TT|) and to two models which has a resemblance to the Cu02 
plane of the high temperature superconductors (chapter VIII.). In the latter case two sim- 



plified models are presented where the solution of vertex equations are essentially different. 
Finally, in chapter IX. we give a brief comprehensive conclusion with some hint concerning 
the applicability of the model. The chapters VII. and VIII. can be read independently. 

II. THE MODEL 

The model proposed consists of two electron bands in the tight binding approximation: 
the broad light (/) band with bandwidth D and the narrow heavy {h) one with energy e^ 
{\sh\ ^ D) measured from the Fermi energy £p. The /-orbitals are on each atom. The heavy 
orbitals are located at some of the atoms in the cell and they overlap only with /-orbitals of 
the nearest atoms. In general theory they may be more than one heavy orbital at one site, 
which corresponds to index 7. In Fig. || we are showing such a model in one dimensions, 
where the light orbitals are s orbitals and and the heavy orbitals are of p and d type. In 
Fig. 1^ as an example of a two dimensional model, the CUO2 plane and the apical oxygens 
below the Cu are shown. The heavy orbitals are associated with the p^ and Py orbitals of 
the apical oxygen. 

The hopping Hamiltonian can be given in the terms of the annihilation operators Cna^, 
and hn-yai where n stands for the position of the atom with heavy orbitals, the spin is cr = ±1 
. The (n, 6) labels those light atoms, which are the first neighbors of the site n at (n + 5/2)a, 
or at site n for 5 = 0, where a is the lattice constant. The labels of the c-operators are not 
defined in a unique way. For example in the case of square lattice both (n, 8) and (n+5, —5) 
labels the same atom on the positions (n + 5/2)a. Thus the one-particle Hamiltonian Hq is 

n,(T n,cr,7 

+ 5Z tlshU.c^Sa+ Yl WcLc„5V + h.c. , (2.1) 

where es is an energy splitting, tgs' and tj^g are hopping parameters (^^^=0 *^ ^ss'^ ^^^ 
the definition n^g^ = c'^g^c^s^ is used, t^ g must hold because of the different symmetry of 
the orbitals. The hopping tgg' in the light band may include the site 5 = as well. The 



direct weak hoppings between the light and heavy bands will be taken into account as a 
broadening of the /i-orbitals (F), which may serve as a low energy cutoff in the logarithmic 
integrals. The part of Eq. ( |2.1| ) due to the /-orbitals can be diagonalized and only the band 
crossing the Fermi energy is kept. In the new band the annihilation operator is denoted by 
(iko- where k is the momentum. The contributions to the Fourier transforms of the Cka^ of 
the d-band crossing the Fermi surface are 

CkSa = 06(k)cik<7 , (2.2) 

where 05 (k) are the amplitudes of the electron on the orbital denoted by d in the state k 
and the contributions of the other bands are dropped. 

The interaction Hamiltonian consists of two parts, H[^t = Hu + Hi. The Hamiltonian 
describing the Coulomb repulsion is given by 

^^ = E E U2hl^J^,^,c^,^,h^^^ + h.c. , (2.3) 

n,6 7,o',(t' 

where U^ is the Coulomb integral between the heavy orbital 7 and light orbitals at site 5 
(see e.g. Fig. |l]c). The Hamiltonian due to assisted hopping is 

n,5 7,cr 

= E E iyisAs-.^^r^s-.K^a + h.c. , (2.4) 

n,<5 7,0- 

where t] is the amplitude of the assisted hopping between the heavy orbital 7 and light 
orbitals at site 6. The possible role of the on-site correlations between the heavy electrons 
are taken into account in the renormalized parameter e^ which can be used if only the single 
excitation of smallest energy is considered at each /i-site and the other excitations shifted in 
energy due to the large on-site Coulomb interaction are dropped. 

The form factor responsible for the orbital Kondo-effect arises from the Fourier transform 
of the Hamiltonian Hi describing the assisted hopping, 

Hi = w2 E E^>"'^'^^''"''^''^^''^^(kO0Kk2)0,(k3) 

:<X.-A3->7k.<x + h.c. , (2.5) 
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Ar2 

ki,k2,k3 5,7,0- 



X( 



where k/i = ki + k2 — ks, the notation given by Eq. ( |2.2| ) is introduced and A^ is the number 
of unit cells in the sample. This form factors prevents the cancellation of the logarithmically 
divergent loops which usually cancel if the interaction terms are diagonal. In Fig. ^ the 
diagrammatic representation of the bare interaction is shown. 

Along the calculation several approximations will be applied. The momentum integrals 
are split to an energy integral perpendicular to the Fermi surface and a momentum integral 
on the Fermi surface. Thus e.g. in d-dimensions 

where dS^, is a surface element of the Fermi surface, k is the unit vector directed from the 
center of the Brillouin zone to the surface element dS, and k^(K;) is the Fermi wave vector in 
the direction k. Furthermore, the form factor appearing in Eq. ( p.5|) and the wave functions 
05 (k) are replaced by their values taken at the Fermi surface like 

g-i(ki + ks - k3)<5a/2 _^ ^-i[kp{K,i) + kp{K,2) - kp{Ks)]6a/2 / ^^2.7) 



and 



(k) ^ <Psm^=.(.^ = M>^)- (2.8) 



lk=k^(K) 



The Fermi velocity Vp{K,) is defined as 

de 



dk 



= Vp{K) . (2.9) 

k=kp(K,) 



III. VERTEX CORRECTIONS 

The renormalization scheme consists of two steps due to: 

(i) Coulomb interaction Hu which in the leading logarithmic approximation results in 
the corrections tU^ln"' \D/e\, (n=l,2,3...) to the vertex t (see Fig. ^), where the energy 
cutoff e is the largest of the energy variables, e = Max(|u;|, \eh\, F); 
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(ii) assisted hopping Hi in the next to leading logarithmic approximation by calcu- 
lating the self-energies proportional to ut^ In D/e and the /-particle four vertex function 
~ t^lnD/e. As a consequence an identity of Ward type, the self-energy depicted on Fig. |^ 
is connected with the vertex correction giving the first contribution to the induced interac- 
tion V in the electron-hole channel (see Sec. V). 

In both steps the /i-orbitals at different sites contribute independently. Joint contribu- 
tions of different sites are relevant only in higher logarithmic orders which are neglected. 
The intermediate states with two /-particles or two holes do not contribute to the renor- 
malization as logarithmic contribution arises from the particle-particle channel only if their 
total momentum is zero. The latter is the case of the superconducting gap equation, but 
not in the processes described above. 

The strong renormalization of t is generated by three diagrams which are depicted in 
Fig. P?. The summation of these diagrams are performed in the leading logarithmic ap- 
proximation considering the "parquet" diagrams. As t is a small variable compared with U 
(|t| ^ |/7|), the equation are linearized in t; thus two or more /i-particle intermediate states 
are not included. This approximation is consistent as the parquet scheme does not contain 
these diagrams. The parquet equation for vertex t can be written in a schematic form as 

r^ . 1 

t{uo) = t + --- t{e)- — -U{e)de, (3.1) 

JMax(|a;|,|£fc|) [-^) 

where the momentum integrals and the coefficient are not presented and U{uj) = U must be 
taken as it will be discussed later. The variable e represents the smallest energy denominator 
in the one /i-particle and one c?-particle (hole) channel. 

The differential form of the schematic Eq. ( |3.1|) is the scaling equation 

J-^ = ---i{u)U. (3.2) 

OUJ 

The renormalization procedure generates an effective assisted hopping, which can be 
given as 

7,0- ki,k2,k3 
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where k/i = ki + k2 — ka. The unrenormahzed hopping [see Eq. ( p.5|) ] is 

s 
Similarly, the niomentuni representation of the effective Coulomb repulsion is 

■JpY. Yl U\K.^,l^2)hi^,^A,-ad^.,-aK,,^a- (3-5) 

7,cr ki,k2,khi 

where k/,2 = k/^i + ki - k2. 

The present treatment is essentially simplified by the fact, that logarithmic vertex cor- 
rection of the Coulomb vertices shown in Fig. ^ cancel out as they describe static potential 
represented by the /i-particle, thus the scattering strength Uj is not renormalized. Thus the 
effective Coulomb interaction U'^{n, n') remains diagonal and can be given as 

f/^(Av, k') = ^ f/Je-*['^^(")-'^^("')]'''^/V5(^)05(A^') . (3.6) 

5 

In the case of the assisted hopping such cancelations does not occur, since the form 
factors in front of the loop contributions are not the same, and we can get the following 
scaling equation in leading logarithmic order: 
d ~ 

a;— t^(Kl, K2, Ats, ^) = p{U^{t^l, K)t'^(K, K2, K3, Uj)) ^ 

+ p{U^{k2, K)r{Kl, K, K3, ^))^ 
- P{i'^il^l, 1^2, 1^, Uj)U'^{k, Kg))^ , (3.7) 

where the three terms on the right hand side corresponds to diagrams shown in Fig. ^b. 
Here we introduced the notation of the average of a function /(k) over the Fermi surface as 

ifi'-)).= 77^ f^fi'^)^ (3-8) 



{27tYp 



Vp[K,) 



'^WrJ^- <'■'' 



in dimension d, and p is the density of states, 

Qj I CltJ 1^ 

M 

The scaling Eq . ( |3.7| ) can be transformed to a real space representation by assuming the 
Ansatz 
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i^{Ki,K2,K3,Uj)= Y^ il^g^s^{ij)4>*s^{Ki)^*s^{K2)4>5.^{H,3) 



^1,62,63 



xe 



-i[lip{K.i)Si+\ip{K2)S2-k„{K.s)S3]a/2 



(3.10) 



(where 6 points to nearest neighbor or (5 = 0), thus the following closed set of equations is 
obtained: 



iU- 



S1S2S2 



(uj 



duj 



P 



Y.{KFsJls2S.,i^) + UlFs,sn,,,,^{uj) - ils2si^)Fss,Ul) . (3.11) 



Here the Fggi incorporates the form factors by the definition 



Fss' = (e^*^- 



ikp{K.){S~S')a/2 



(ps{l^)(^*S'{l^)) 



(3.12) 



which appears in the transition amplitude of an /-electron between two sites, 6 and d' , in 
the tight binding approximation. As a special case, the ps = pFss gives the partial density 
of states at a site 6. 

The solution of the differential Eq. ( p.llD can be given immediately as: 



ts^S2Sii^) 



E 



UJ J 



SiS{ 



iT' 



f7(0) 
'^ xf xf sf 



(3.13) 



S'-,S: 



s^s 



where we used the notation F and U''' for the matrix F^g' and the diagonal matrix Uj^, = 
UgSss'- Unfortunately this form of the solution is not suitable for further calculations. In 
the actual computation let us consider the eigenvalue problem of the matrix U'^F appearing 
in the exponent: 



Y.UlFss's]f = X]s 

S' 



7J 
S ' 



(3.14) 



where j labels the eigenvalues thus the spectral decomposition of the matrix U'^F, which is 
not necessarily uniter, can be given as 



THP , — \^ \7„7J»,7J 



(3.15) 



with the row vectors r orthogonal to the column vectors s, thus 
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Y^r^' s-^' = 5,, . (3.16) 

In general the matrix U^F is not invertable, that is why the left and the right eigenvectors 
can be different, such an example is presented in Sec. VIII. (case A). 

Using the identity (U^F)^^/ = (FU^)*,^, which follows form Eq. (|]T|) and that U^ is 
diagonal and real, thus the left eigenvectors of FU^ are the s"''^* and writing tl § § as a 
linear combination of the eigenvectors s as 

n.s.s. = Y.%^)€^'^" ^ (3-17) 

ijk 

the scaling Eq. ( p.ll| ) can be diagonalized and we obtain a set of decoupled linear differential 
equations for the tj^k^oj) 

^^7^ = (A. + A, - \k)p%,{u^) , (3.18) 

whose solution is 

*;-=(-) =€'(^) ■ (319) 



Putting Eqs. ( p.lO| ), ( |3.17| ) and ( p.l9| ) together, we get the fully factorized form of the 
effective assisted hopping: 

t^(A.i,K2,AC3,0;) = $^C(^l)e;('^2)CM(§) '^^ '"tj?, (3.20) 

ijk ^ ^ 

where we have introduced 

eJ(A.) = Y. 4>l{K)sye-^'^^^^^'^/' . (3.21) 

s 

The exponents are of the order of pU, since in the eigenvalue Eq. ( ^.14[ ) F is of the order of 

unity, thus the magnitude of A is determined by U. 

The couplings f^-^ can be determined from the initial (unrenormalized) couplings (with 

D being the bandwidth) i'^^^^Ki, K2, ks) = i^i^i, ^2, i^s, oj = D) [see Eqs. ( |3.4|) and (|3.2CI| )], 

since in that case 
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i^^'\K„K„K„D) = j]c(«i)e;MC('^3)t;?^ . (3.22) 

ijk 

After some algebraic manipulations we get 

ilS^ = J2^lrrrrr7''- (3.23) 

s 

The latter equation tells us that i'jjj is symmetrical in the first two indices ( ijjj = ij^jj) 

and as a consequence [see Eq. (|3.2CI|) ] the t'''(Ki, H2-, 1^3,1 ^) is also symmetrical in the first two 



K variables: 

t'^(Ki, /ta, «;3, w) = ^'^(^2, «;i, «;3, w) . (3.24) 

For later use we mention, that the average of ^*^ can be expressed as 

A similar approach was used by DzyaloshinskiS in case of the two dimensional nearly 



nested electron gas. However, unlike in our case, those scaling equations do not form a 
closed set. 

IV. LIGHT PARTICLE SELF-ENERGY 

Unlike to usual renormalization procedure, where the self-energy corrections become 
important only in the next to leading logarithmic order, in our case we expect a large mass 
renormalization for the light electrons as the logarithmic corrections in the self-energy and 
in the generated electron-electron interaction occur on the same level. 

The /-electron self-energy diagrams are shown in Fig. |^. We are going to evaluate them 
for those electrons which are close to the Fermi level and we assume that \eh\ ^^ sp or for 
energies close to Eh- The latter case is of no importance for later calculation and we are 
presenting it just for being interesting. 

The heavy particle self-energy of the lowest order contains three /-electron lines in the 
intermediate state, thus it is nonlogarithmic except in one dimension. 
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A. Self— energy near the Fermi-level 

In the logarithmic approximation the intermediate state with the two smallest /-electron 
energy variables ^i and ^2 must be picked up and then both vertices must be fully renor- 
malized, but their lower cutoff must be replaced by \eh\- A typical integral to be calculated 
for the self-energy with energy 00 <^ Eh is 

Je, ^'Js, ^'Ui + 6; ±sign(e,V-(ei + 6 + k.|) 

= -{u;Teh)A{D/\eh\), (4.1) 

where the sign is appropriately chosen to give the correct contributions of different diagrams: 
the upper sign stands for the first two diagrams and lower sign for the last in Fig. §. In 
the actual calculation we pick up the most divergent term, so that the largest exponent 
is taken as a. The vertex contributions contains also u variables, but its sign is not well 
defined. It is easy to demonstrate that the appearance of such an u changes the right hand 
side of Eq. ( |4.1| ), by factors like (1 ± a)~^. As the logarithmic approximation is applied, 
therefore, the corrections proportional to U must be neglected. Furthermore, the shift of 
chemical potential EhA is also only an approximate value, since the in this approximation it 
is difficult to determine exactly^. The coefficient A can be calculated and 



A{D/\s,\) = - 
a 



D 



(4.2) 



is obtained. 

The complete expression for the self-energy is 

-D 

'Max(<^,£h) 



I'D 

JMa,x(uj.eh) 



- (t\K,\n\n;e)-t^*{n\n'\K,-e)) J , (4.3) 



16 



where the different terms on the right hand side comes from the three diagrams in Fig. 
respectively. Performing the integration, we get 



+ £h[xi{i^) + X2(k) - X3{i^)]p'^A{D/eh) , 



(4.4) 



where 



X2(k) 

X3(«^) 



E 






(4.5) 



and, furthermore, the most divergent term of P [see Eq. ( 3.20 )] is used (the exponent a is 
the largest of the eigenvalue combination Aj + Xj — A^), which is a good approximation for 
most of the cases. 

The symmetry of P [see Eq. ( |3.24 )] implies Xil**) = X2{i^), furthermore, the average 
over the Fermi surface of every Xji^^) is equal and will be denoted by x- 

Knowing the self-energy given by Eq. ( [4.4|) , the renormalized one particle Green's func- 
tion is obtained as 



G,(k,uj) 



(4.6) 



^-^F,rcn('«)[^-^F('«)] 

where we have neglected the renormalization of the chemical potential, which comes from 
the second term in Eq. ( [4.4|) and it is hard to calculate in leading logarithmic order approx- 
imation. The renormalization constant is defined as 



^-1 



1- 



dReEi(ij,k) 



duj 



u)=0 



l + [2xi(/^)+X3(^^)]p'^(^/k^|) > 1 



and the Fermi velocity is renormalized nearby the Fermi surface as 



(4.7) 



F,rcn 



K 



Vp{n)Z^ < Vp{n) , 



(4.8) 
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thus it is suppressed, leading to mass enhancement. The dispersion curve is schematically 
plotted in Fig. |^. In the case of the large renormalization described by Z^ ^ 1 the large 
modification of the dispersion curve is expected for |ci;| < \eh\ and the renormalization 
gradually disappears as approaching larger energies, |ti;| > \eh\- Thus, the large enhancement 
of the density of states is restricted to a small energy region \uj\ < \eh\Z^. 
The average strength of the mass enhancement can be given by 

(Z,i) = l + 3x/A(D/|e,|). (4.9) 

The mass renormalization can be very large, thus the scaling may result in heavy 
fermionic behavior. The calculation presented are justified only in the region a <^ 1, but 
there is no indication for nonanalytical behavior, thus in order to get qualitative result the 
extrapolation for intermediate coupling a < 1/2 is adequate. As it is known from the X-ray 
absorption problemQ, for large Up that quantity must be replaced by 5/it where S is the 
phase shift {6 < 7r/2). 

In Eq. ( [4. 61 ) there is also a shift of the energy proportional to Xil**) + X2('^) — Xsii^) 
which is sensitive on the direction k. That shift changes the number of the electrons inside 
the Fermi surface, thus the Fermi energy must be corrected in order of p'^x^h^iD/ehjZ,, to 
keep the number of particles inside the Fermi surface unchanged and that is associated also 
with the deformation of shape of the Fermi surface. 

Finally it must be mentioned, that even a large wave function renormalization does not 
play an important role in the vertex equation for t, as the Green's functions are taken 
between two points in the real space not further apart than the lattice constant. In that 
case the Fourier transform of the Green's function multiplied by a slowly varying function 
like exp(zfc^.a/2) must be integrated with respect to the momentum. The result is almost 
independent of the self-energy and the correction is 0{U / D). 



B. Self— energy near the singularity 

If the energy of the electron is large enough, than the scattering of an electron from the 
filled Fermi sea to the unoccupied heavy level becomes a real process, and we are faced with 
a problem very similar to that of the X-ray edgeQ, where a power-law behavior is observed 
in the absorption function. In our case the role of the X-ray is taken over by the incoming 
light electron, and power-like behavior is expected in the imaginary part of the self-energy 
and so in the spectral functions. For simplicity we are going to consider the Eh > case 
here. Similar considerations are valid if Eh < 0. 

It is easy to see, that in Fig. || for u; ~ Eh the first diagram S/^i(c<j, k) and second diagram 
S;^2(^! '*); while for u ^ —Eh the third diagram £^^3(0;, k) is singular. In this case the integral 
in §3) is 

J\.-e,\ "^^^ J\.^e,\ "^^^ Ui+6; ±^-(^l+6+kh|) ^ 

= -{uj T eh)A{D/\uj T eh\)[l ± ia7iQ{±uJEh)] , (4.10) 

where the G(x) is the step function and it ensures the proper analytical behavior. The signs 
are the same as they were in Eq. ( |4.1| ). 

For this special choice of energies the renormalization of the assisted hopping [Eq. ( p.?! )] 
should be considered more carefully. It turns out that the energies of the internal electron 
lines are such, that for S;^i(ci;, k) the first term in Eq. (|3.11|) is not singular, similarly for 



S;2(^,'*) the second term and for the T,i^3{uj,k,) the third term can be neglected, so that 
the exponent Aj + \j — Xk is reduced to the values \j — Xk, Aj — A^ and Aj + Xj respectively. 
The exponent a should be associated with the largest exponent for each case separately and 
will be denoted by ai = 0:2 and a^, and Ai is the function defined by Eq. ( |4.2| ) with a = ai. 
Putting everything together, the self-energy for energies near cj = e/j is 

S,(k; uo) = -{uj - Eh)[xiii^) + X2{t^)]p^Ai{D/\uj - Eh\)[l + iai7TQ{u - Eh)] (4.11) 

and 
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S^A^; uj) = -{uj + eh)xz{i^)p'A^{D/\uj + eh\)[l - la^TiQ^-oj - sh)] (4.12) 



for uo K, —Eh- The x(k)-s are the ones defined by Eq. (|3.21| ), however the terms to be 



averaged are not the same and, therefore, the averages are no more equal. 



The imaginary parts of the self-energies (|4.11|) contribute to the spectral densities of 



the particles and it emerges as soon as the energy of the electron is larger then ^ Eh-, while 
the self-energy ( [4.12| ) gives similar contribution for the spectral density of the holes having 
energy below the threshold —Eh- 

V. GENERATED INTERACTION BETWEEN THE LIGHT ELECTRONS 

To address the possibility of superconductivity in light band in these models, we need 
to calculate the generated interaction in this band for parallel- and anti-parallel spin. The 
typical diagrams are shown in Fig. |^. In the logarithmic approximation the intermediate 
state with one /i-electron and the /-electron with the smallest energy must be picked up 
and in vertices that energy must must be used as the lower cutoff. 

The schematic form of the generated interaction is 

t{E) t{E)dE ~ -t^{E) (5.1) 

without numerical factors containing the strength of the Coulomb interaction. 

The calculation is straightforward in the time-ordered diagram technique. For the inter- 
action in the anti-parallel spin channel the contribution of the four diagrams in Fig. |^(a) 
is 

-D 



V^{Ki, K2, K3, 1^4) =-p^ 



dE 



^ - Ma.x{ui,eh) 
-{i^{Ki, K, K3; E)-i'^*{n4, K, H2; e))^ - {i^{K,2, K, K^, E)-f'*{n^, K, Hi] e)) \ . (5.2) 

The evaluation of diagrams shown in Fig. ^(b) gives in the parallel spin channel 
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-D 

de 



V^^Ki, K,2, K3, K4) = -T^y] / 



e e 

-(t'^(K, Ki, K3; e)-i'^*{K, Hi, K2] £))^ - {i'^{i^, 1^2, 1^4] e)-i'^*{n, Kg, Ki; e)) J , (5.3) 

where we have antisymmetrized the interaction. These expressions can be evaluated by 
using Eq. (p.2CI|) for the scattering amphtude. The induced interaction is constant in the 
energy range |ci;| < Z\eh\ and drops rapidly for higher energies. 

The interaction can be split into two terms: the interaction V^ between electrons forming 
a singlet and V'^ between electrons in a triplet state. In the singlet state the spin part of 
the wave function is antisymmetrized, so the V^ is symmetrical in the ki and H2 variables 

V^{ni, K,2, Ks, 1^4) = V^{k2, Ki, K3, K4) . (5.4) 

For the triplet case the symmetrycal spin wavefunction implies antisymmetrycal property of 
the V'^, so 

V'^{Ki, K2, K3, K4) = -V^{k2, Ki, K3, Ki) , (5.5) 

and similar holds for the last two k, variables. Furthermore, it is evident that for both 
interactions 

V'^'^iKi, K,2, Ks, ^^4) = ^^'^{^,2, Kl, K4, 1^3) ■ (5.6) 

The V^^ is by definition the interaction between electrons in triplet state. To V^, however, 
both the singlet and triplet states contributes. We can split them by using the symmetry of 
the V^ and V'^: 

^•^(Ki, K2, K,3, Ki) = -[V-^{Ki, K2, K3, K4) + ^"^(«^2, «^1, 1^3, 1^4)] , (5.7) 

and 

V'^{Ki, K2, K3, K4) = -[V-^{ki, K2, K3, 1^4) - V-^{k2, Ki, K3, K4)] . (5.8) 
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The latter must be exactly V^^{k,i, k,2, /^a, ^4), as it can be seen if we take into account the 
symmetries of f^ [Eq. (|3.24|) ]. This separation is possible because of the absence of the 
spin-orbit coupling in this model. 

Performing the energy integral, we get for the interaction in the singlet channel 

7 
-{f'{ni, K, K3)t'^*(K4, K, K2))^ - {^{^2, K, K4)^'^*(k3, K, l^l)) ^ 

-ip{n2, 1^, l^3)i^*{K^, K, Ki))^ - (^(Ki, K, K,4^)i^*{K,3, K, ^2)) J (5.9) 

and for the triplet channel the result is identical with Eq. (|5.3| ) where the energy integral 
must be replaced by A{D/eh)- 

We can see that if the assisted hopping i'^ is structureless, than the generated interaction 
disappears both in singlet and triplet channel. However, a small k dependence can lead to 
finite interactions which may be enhanced due to large A{D/\eh\)- 

VI. SUPERCONDUCTING TRANSITION TEMPERATURE 

As it has been mentioned in Sec. II the induced interaction between the light particles 
(Fig. P) and the self-energy (Fig. H) are connected by identities similar to the Ward identities. 
The Ward identities are valid only in the electron-hole channel with zero total spin and not 
in the Cooper channel. These relations are due to the fact, that both occurs first in the 
second order of the perturbation theory on t. That is in contrast to most of the other 
logarithmic theories, where the vertex correction exists already in the first order. Thus in 
the present case, the induced vertex and the self-energy must be treated simultaneously in 
determinations of the superconducting order parameter A and of the transition temperature, 
which are based on the selfconsistent equation shown by diagrams in Fig. ||. 

The order parameter is defined as Ao-o-'(k) = (ck,CTC_k,o-') , where a and a' stand for spins. 
Separating the spin wave function, A is either even or odd function of k, depending whether 
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we are looking at singlet (A'^) or triplet (A'^) Cooper-pairs. Using the interaction V^^'^\ 
the self-consistent equation for A'^^-'") is 



A^(^)(k) = -T.l Yl E ^"^^^k, -k, -k', k') 

k' U!,n 

X G(k>„)G'(-k>„-^™)A^(^)(k') , (6.1) 



with Tc being the transition temperature and Um = (2m + 1)ttT are the Matsubara frequen- 
cies. After summing over the internal energy, we arrive at 



A^(^)(..) = -T^p J2 A(^"^^H-> --> --'> -')^(-')A^(^)(..')).^ , (6.2) 



UJm<l^. 



m^^c 



where we introduced the frequency cutoff Uc as we have replaced the Um dependent quantities 
by their low frequency values, assuming that Tc is much smaller than the characteristic 
energy Eh- In our case the Uc is determined by the fact, that the contributions for the 
integral over the internal energy comes from energies smaller than Z\eh\, as it is determined 
by the linear part of the dispersion curve where the density of states is enhanced (Fig. ^, 
so the summation over the frequencies yields 

A^(^)(A.) = -In (i^^^) p{V'^^\>^,->^,->^',>^')Zi^')A'^^\^'))^, . (6.3) 

In other words, in the selfconsistent equation for the order parameter each interaction 
V is associated with two Greens's function. Only the contribution of electrons with energy 
uj < Eh is kept as the vertex drops rapidly for larger energy. For the dominating low energy 
electron the strength of the poles is Z^., and its energy \uj\ < ehZ^ (see Fig. 0). Furthermore, 
the density of states p for such electrons is also enhanced by Z^^ [see Eq. ([4.8|)] . 

In the singlet channel, where it can be assumed that the K-dependence of the super- 
conducting order parameter A(k) is weak, the approximate strength of the dimensionless 
effective coupling is obtained by evaluating Eq. (|67^ ) so that 



gip r^ p{V'{K,-^,~^',^'))^^^,{Z{^'))^, 

r^{Z-Y'{V')p, (6.4) 
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where in the last approximation Z^, is replaced by its average over the Fermi surface [see 
Eq. ( |4.9| )]. If the K-dependence of the V and Z is not large, then the average of V^ is of 
the order of px^- Inserting the calculated value of Z^ given by Eq. (^^91), we can see that 
the dimensionless coupling is 

where q is of the order of unity and depends on the details on the /^-dependencies, which 
is approximated in Eq. (|6.4|) . If the mass renormalization is large, 1 + 3xp'^A{D/6h) ^ 1, 



(S) 

then the denominator is important and g^f^ saturates. 

In this way in the BCS theory the transition temperature is 

Te = |£,|(Z,)e-i/^^ff' (6.6) 

which can be expressed as a function of the averaged mass enhancement {Z~^) given by 
Eq. ( [4.9[ ) and of the parameter q. That function is shown in Fig. ^. The transition tem- 
perature Tc is suppressed by the decreasing (Z^) for strong renormalization. The critical 
temperature is the largest in an intermediate region of moderate renormalization where the 
mass enhancement (Z^^) ~ 2 — 5. As we have seen, in general it is essential to include the 
renormalization of the Fermi- velocity due to self-energy to get a consistent treatment of the 
superconductivity. 

VII. THE ID MODEL 

To demonstrate the procedure described in the previous sections, let us apply it to a 
1-dimensional case introduced by one of the authorsO. In this model the light orbitals are 
on every site n and the heavy orbitals are on every second site, as shown in Fig. |1|. The 
degeneracy of the heavy orbitals is neglected thus the index 7 will be dropped. The hopping 
to heavy orbitals is t^ ^ = th and the hoppings between the neighboring light orbitals are t, 
the Coulomb repulsion between the light and heavy orbitals is 
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f/ , a 6 = 

Us = { (7.1) 

u' , a 6 = ±i , 

where the on-site and nearest neighbor interaction is kept. Finally, the assisted hopping is 

f , a 6 = 

ts={ (7.2) 

It, a 6 = ±i . 

in case of light s and heavy d orbitals, where the wave functions have the same sign in the 
regions of the overlap. In the previous workllll U' = was taken. 



The dispersion relation of the fermions in the light band, shown in Fig. 10, is given by 



Ae / / Ae: \ , „ ^ ka 



uj^(k) = — ±M—j +4t2cos2- (7.3) 

in the tight binding approximation, where the +(— ) stands for upper(lower) band, the 
Brillouin zone is — vr/a < k < ir/a, and Ae = es — eQ is the relative shift of the atomic levels. 
As a next step, we have to determine the amplitude 0(k) for 6 = ±1,0 [see Eq. ( |2.2| )]. 
In the case of even tss' the 0-s are not k dependent and 0+i = 0_i. Since the probability 
is normalized, 0^^^ + (p^ = 1, we can parametrize (p-s by 0+1 = 0_i = cosip and 0o = simp. 
Here the subscript +1 (—1) stands for the right (left) neighbors for the heavy orbital. If the 
energy levels of the light orbitals are the same (es = 0), then all the 0-s are equal to l/-\/2 
[if = IT /A). In the general case the ip can be estimated as 



Ae IfAe^^ 



tan(^ = ^ + W ^ +1, (7-4) 

^ 2tcos?9 y V2tcos^y ^ ' 

where d = kpajl is introduced in accordance with Ref. |1T[ It depends on the filling i/, 

"d = 11^/2 for the lower band (z/ < 1) and d = 7r(l — z//2) for the upper band (z/ > 1). 

Specially, in one dimension, Eq. ( |2.6D reduces to 



where k = 1 and —1 denotes right and left moving electrons, with kp{l) = kp and kp{—l) = 
—kp, similarly for the velocities Vp{±l) = ±Vp. The form factor appearing in Eq. ( f^.7D is 
also simplified significantly, since there are only four combinations exp(±z7?) and exp(±z3'(9). 

25 



The density of states p [Eqs. ( |3.9| )] is the usual one for the tight binding approximation: 

1 2a 



P 



2tt vp 



(7.6) 



and the calculation of the amplitudes F^s' [Eq. ( |3.12|) ] gives 

Foo = sin^ ip 

F++ = F = cos^ y) 

Fo+ = Fo_ = cos "d cos ip sin p 
F+_ = cos 2'd cos^ p) . 

The UF matrix is then 
/ 



(7.7) 



UF 



U' cos^ p U' cos d cos v? sin p U' cos 2'i9 cos^ p 

U cos -(9 cos p sin (y9 [/ sin^ p U cos "(9 cos p sin </) 

, t/' cos 2'd cos^ y) f/' cos d cos v? sin p U' cos^ (/) , 



(7.8) 



The matrix UF is not symmetric, thus the left and right eigenvectors are different, and one 
obtains for right and left eigenvectors 



=(1) 



U' cos p cos t9 

U sinp 
U' cos p cont'd 



=(2) 



smyj 

—2 cospcos'd 

smp 



J3) 



1 




(7.9) 



and 



-(1) 



-(2) 



-(3) 



1 

1 



'cosy^cos-i? , sinp , cos </) cos -(9) 
{Usinp, — 2f/'cosv9cos'i9 , Usinp) 



2Ai 

(1/2, 0, -1/2) 



(7.10) 



with eigenvalues Xi = U sin^ p + 2U' cos^ p cos^ -(9, A2 = and A3 = 2U' cos^ y? sin^ {}. 



The functions ^ [see Eq. ( p.21|) ] are 



6(/t) = Ai 

Uk) = 

,^3(fi;) = — ^2 sin K-i? cos (y? 



(7.11) 



26 



The nonvanishing averages are 





iU^) ') = A? 




(|^3(ft:)2) = 4cosVsin'^. 


Eq. d^), we get 




f(0) 

''111 


= 2t--3 cos^ 'd cos^ (f 


f(0) 
''133 


r(0) 7(0) t 1 
= %3 = 431 = 2;^ COS^COS(/? 



(7.12) 



(7.13) 



so that for the effective assisted hopping we get 



t{Ki, K2, %; a;) = 2t cos^ ipl [cos^ d + cos d{sm Kid + sin K2'd) sin K^^-d] {D/ujY ^ 

-cosi9sinKi^sinM(^/c<^)^^^^'"^'^} • (7.14) 

Since in the assisted hopping only the transitions of an electron from one of the neighbor- 
ing sites to the heavy orbital play role, the amplitude cos^ y? is easy to understand being 
associated with the three light electron lines. Furthermore, we can see that the Coulomb 
interactions appear in the exponent only as the combinations U' cos^ y? and U sin^ 99, what 
means that the Coulomb repulsions enters only as effective repulsion normalized by the 
single site fermion densities. 

As K can have the values ±1, the assisted hopping can be given by the four amplitudes 
ti, ^2, ^3 and ^4, introduced in Ref. |n|. They can be expressed by t(fi;i, K2^ k^; uj) as 

ti=t(l,-l,l;^)=r (-1,1,-1;^) 
t2 = t(l,-l,-l;^) = t*(-l,l,l;w) 
t3 = t(l, 1, -1; 07) =r (-1,-1,1;^) 
ii = i{l,l,l;uj) =r (-1,-1,-1;^) (7.15) 

so that 
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ii =^2 = 2tcos'^ V9 
ii = 2icos^ if 



cos^ ^ - 2 cos ^ sin^ ^) {D/uY^' - cos ^ sin^ ^{D/)"^^^'-^''^ 
cos^ ^ + 2 cos ^ sin^ ^) {D/luY^' - cos ^ sin^ i9{D/ujy^^^'''^'^ 



(7.16) 



These equation are exactly the solutions of the Eq. (7) in Ref. jTll with initial couplings 
^(0) ^ ^io) ^ ^io) ^ ^~cosVcos?9 , 4°^ = tcosVcos3?9, ?7' = and (^ = 7r/4. 

As a next step we calculate the generated interactions between the fermions in the light 
band. We are using the notations common in the theory of the one-dimensional Fermi gas. 
For a review, we refer to the paper of Solyomtl. 

In the antiparallel channel the forward scattering, denoted by g2, is 



r72 = ^^(l,-l,-l,l)=2|tip-2|t3 



(7.17) 



Replacing the assisted hopping by its most divergent part, determined by the largest expo- 
nent, we get 



3 ,„\2 



g2 = (2t cos (f) 



2sm.'2dcos2'd{D/ehy^^ , if Ai > A 



(7.18) 



, if Ai < A3 . 

If Ai > A3, than (72 is positive for small fillings and became negative for larger fillings. 
For the backward scattering, where the momentum transfer is large {2kF), we get 



gi± = v^{i, -1, 1, -1) = 2ti(t2 - h) + 2ii{i; - i 



4; 5 



(7i|| = V\\l, -1, 1, -1) - V\l, -1, -1, 1) + ^2 = 4|tip - 2i% - 2UI . 
Since ti = ^2) the gi_[_ = gi\\ = gi holds, which means that the model is isotropic, and 



(7.19) 



3 ,„n2 



gi = (2tcos if) 



sin^ 2^(1 + cos 2^){D /Eh 



,2Ai 



if Ai > A3 



sin2 2?9(l - cos2^)(D/£;,)4^3-2Ai ^ if Ai < A3 



(7.20) 



They are presented in Fig. |TI. 

In one dimension the values of these couplings determine the nature of the ground statea. 
We analyse the Ai > A3 and the A3 > Ai region separately. In the first case the gi is 
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negative, and depending whether the sign of the combination gi — 2g2 is positive or negative, 
the expected ground state is singlet superconductivity or a CDW. A simple calculation 
shows that for § < dc the ground state is CDW and for larger filling {d > dc) we expect 
superconductivity, where cos2dc = —1/5 (this corresponds to filling nc = 0.564). On the 
other hand, if the exponent 2A3 — Ai is larger, then gi — 2g2 is always positive. Since gi 
is positive, we get triplet superconductivity. The two exponents are equal if U sin^ ip = 
—2U' cos^ (f cos2'i9. The phase diagram for this model is then shown in Fig. |l2|(a), where we 
can see that the nature of superconductivity changes as the Coulomb interaction increases 
between the heavy orbital and light orbital on the neighboring site. 

There is also a strong mass renormalization in this model. It is easy to see, that xi = 
X2 = Xs = X ^'I'e equal [see Eq. ( ^4.5| )], and 

X~tl + tl + tl + tl. (7.21) 

Putting in the actual expression for t-s, we get 

. , f 4cos2'(9(cos^^ + 2sin^'^) , if Ai > A3 
X = {2tcos^^f < (7.22) 

I sin^2^cos2^) , if Ai < A3 . 

which become large for intermediate fillings. 

If the parity of the light and heavy orbitals is not the same (e.g. light s— or d— orbitals 
and heavy / orbitals, see Fig. 0(b)), then the assisted hopping is odd and 

{0 , a 6 = 
(7.23) 
6i, a 5 = ±l . 

The parity of the assisted hopping has no effect to the eigenvalues and eigenvectors of the 



UF. Using Eq. ( 3.23|) , the t\-l-s are in this case 



7(0) _ i 

'■333 — A 



and for the effective assisted hopping we get 

29 



i^Ss = ifl = if^3 = 4 cos^ ^ cos^ ^ , (7.24) 



t{Ki, K2, K3; uj) = —2itcos^ (p< [sin Ki-i? sin fi;2^sin Ks-i? + cos^ '(9(sin Ki-;? + sin K2'&)]{D/ujY ^ 

-cos^^smK3^{D/uY^^^'~^'^\ . (7.25) 

In this case the ti, ^2, h and t4 are 



pAs 



tl = t2 = -2it cos^ (p sin'' ^{D/luY^' + cos^ 1^ sin ?9(D/cj 



\p{2Ai-A3) 



^3 = —2itcos ip 
ti^ = —2it cos^ v? 



[- sin^ ?9 + 2 cos^ t9 sin ^) (D/cj)^^^ + cos^ ^ sin ^^(D/)''^^^^"^^) 
^sin^ ^ + 2 cos^ ^ sin ^) (D/cj)"^^ - cos^ ^ sin ?9(D/cj)''^^^^"^'^ 



(7.26) 



Comparing them to those of the even assisted hopping [Eqs. [7.16|| , we can see, that disre- 
garding the prefactor of i and the exponents, the only difference is the interchange of sin t? 
and cos'i?. 

Keeping only the most divergent parts in the assisted hopping, for the induced interaction 
we get 



3 ,^\2 



g2 = (2tcos (fi) 



0, 



if Ai > A3 



(7.27) 



-2sin'2^cos2?9(D/£;,)2^i , if Ai < A 



If Ai < A3, then g2 is negative for small fillings and became positive for larger fillings. 
For the backward scattering 



_ ^ ^ , sin2 2^(l + cos2^)(D/e;,)2^3 
gi = {2tcos^^ip) 



if Ai > A3 

sin2 2^(l - cos2^){D/ehY^'-^^'' , if Ai < A3 



(7.28) 



The boundary between the two region with different exponent is the same as it was in 
the case of even is- If Ai > A3, both gi and gi — 2g2 are positive, so the ground state 
is TS. In the region where Ai < A3 holds, the gi is negative, and for -(9 < -(9(7 we get 
triplet superconductivity and for 'd > 'dc the ground state is CDW. Here ^ is defined as 
cos2dc = 1/5 and the corresponding filling is vc = 0.436 . The phase diagram is shown in 

Fig. igb. 
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We can conclude that the main effect of changing the parity leads to exchange of the 
CDW and superconducting ground state. 

In the present calculation we have neglected the interaction between the light electrons. 
That interaction changes the f^-s, so that the phase diagram changes also and new phases 
appear near z/ = and u = 1. 

VIII. TWO-DIMENSIONAL SQUARE LATTICE: Cu02 PLANE 

In the compounds characterized by Cu02 planes a possible representation of the h- 
orbitals are the two non-bonding p-orbitals on the apical oxygens located below or above 
the Cu sites in the Cu02 plane (the 0^~ ions are in the tetrahedral position around Cv^^ 
ion). For details see Ref. p!4| , p!2| . The Cu02 plane forms a two dimensional square lattice 
depicted on Fig. ^, where the /i-orbitals are at the corners and the Z-band is formed by the 
orbitals a^., ay and h at the middle of the sides and at corners, respectively. One possible 
choice of the symmetry of the orbitals corresponds to the Cu02 plane with p^ (flx) , Vy {'^y) 
and dx2_y2 (6), while the two heavy labelled hj •j = x,y are of px- and Py-type . For this 
model the hopping to heavy orbitals is 

tls = t^P's (8-1) 

and the hoppings between light orbitals are 

t..A''" "^'^'"■°' (8.2, 

I t'Pss' 1 otherwise , 

furthermore, 

J.,, if 5' = (0,0) 
65= < (8.3) 

I , otherwise . 

The p^ , p^ g, and Pg stand for the relative signs of the real wave functions in the overlap 
regions. For the nearest neighbor hopping between oxygen and copper in the Cu02 plane 
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Pg = —p_g holds and for the next to the nearest neighbor hopping between oxygens ps^s' = 
±1 with +(— ) sign ii S — 6' is parallel (perpendicular) to the (1,1) direction and for the 
apical oxygen P^ = Pjg if 6 is parallel with 7 = x, y and zero otherwise. 
The Coulomb repulsion is 

{U , ii 6 = 
(8.4) 
U' , otherwise . 

where U is the repulsion between electrons on the apex-0 and Cu, while U' is between 
apex-0 and the next oxygens in the plane. 
The assisted hopping is given by 

n = iP2 ■ (8.5) 

In the following we will consider two limiting cases. 

A. Case of [/' = t' = 

Here we are considering the extreme case where both t' and U' are zero. The dispersion 
relation is given by 

1 



uj± = -[e,±\lEl + AQ^] , (8.6) 

where the +(— ) stands for upper(lower) band, f^^ = t^f, + tyi^ and t^k = 2tsmkaa/2 (a = 
x,y). The 0-s [see Eq. ( p.2| )] are given by 

and (ps = I ., • (8.7) 



The calculations of p and F-s, defined by Eqs. ( |3.9| ) and ( |3.12| )], are straightforward: 



-I^^^¥^'<'W' (-' 



and 
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"00 — VO 



Fxo — Fyo 



f!., -f;,„ 4 (£)(!-) w^ 



1 f 2t\ 

^ XX ^ —x,—x ^ yy yt~y 7) I I l Cj|</)o| 

F =-F =-F =F ='-(^] 

J- xy J- x,~y J- -xy ^ -x,-y O \ , j 



2 r 



K'(c) 



F =F =F =F ='-(^] 

J^x-x J^-x,x ^y,-y ^-y,y \ u I 



2 r 



C - 1 - 2c^ + 2 



l</^0| 

K'fc) 



i0or , 



where c = 1 — i7^/(2t)^ and the elhptic function are defined as 



K{k) 



7r/2 



(i(/9 



V 1 — k^ sm" yj 



"■^ cin , 



and E(A;) = / d(p\Jl — k"^ sm^ ip 



and K'(fc) = K(A;'), E'(A;) = E(A;') where k' = y/l-k^. 
The vectors s and r are [see Eq. (|3.15| )] 



(8.9) 



(8.10) 



=(1) 



= (2) 



—i 




-2/ 



.(3) 






i 
—i 

-V 



.(4) 



.(5) 






1 
1 





f8.111 



and 



,(1) 



,(2) 



,(3) 



.(4) 



,(5) 



{-if , if , -if , if , '^) 
(-V2 , V2 , , , 0) , 
(0,0, -V2 , t/2 , 0) , 
(1/2, 1/2,0,0,0), 
(0,0,1/2, 1/2 , 0) , 



(8.12) 



with eigenvalues Ai = t/|0oP and A2 = A3 = A4 = A5 = 0, and 



/ = (1 - c)t/LU± 



(8.13) 
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furthermore the d indices are ordered as (x, —x,y, —y,0). The ^ functions [see Eq. ( p. 211) ] 
are 






2t 
2—(Posm\kp^{K)a/2) - 2/0o , 

2t 
2—(Posm\kp(K)a/2)-2f^o, 

UJ± 

2t 

— 0osin(/i;^^(K;)a) , 

2t 

— (f)osm{kp (K)a) . 



(8.14) 



The nonvanishing averages of ^$,* functions are from Eq. (|3.25| ) 



{U'^)\') = \M 



{\U'^)\') = {M'^)\') = 2 — 



2t V /E'(c) 



c \(po\ ■ 



(8.15) 



u±J \K'{c) 

Calculating the renormalized assisted hopping, we pick up the most divergent term in 
Eq. (|3.20|) . The largest value of the exponent Aj + Xj — \k corresponds to the choice of 
i = j = 1 and, since the remaining four eigenvalues are degenerate, with k =2,3,4 or 5. 
We can determine the initial values t"^'-"-* from Eq. ( |3.23| ), and we get t^iJ = —ip and 



.2/(0) 
''lis 



-tp with other t^j^ -s being equal to 0. So the assisted hopping is 

2A1-A4 



,x(0) , 



D 






and for P{ki, ^2,^3; 00) the ^4(k;3) is replaced by ,^5(^3), so that 

2f / D\ 'iu\4>a?p 



t^{Ki,K2,K3;^) = -tf 



^''^o|Vo^^sin(A;p, (^3)0 



where hp = {kp^-, kpy) is the Fermi vector parallel to the n. 
The expressions for x-s, according to Eq. (|4.5| ), are 



xi = X2 = t?iV?i4ieip(i6p)(ia(^)p) +iMi5i6p(ieip)(ia(^)i') 
x3(^^)=t?T4tii4(i6r)(i6r)ia(^)r+tMi5(ieir)(ieip)ia(^)r 



(8.16) 



(8.17) 



(8.18) 
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or, inserting the ^-s and tJlJ-s, 



ijk 

E'(c) 



- 1 ^6 ^' 



c' 



P 



1 <?* 

'^ r . 2/7 / N N I • 2/7 / \ mJ2 



_K'(c) 
Xsii^) = ^4^[sm\kp^{K)a) + sin^(A;^^(A.)a)]P . (8.19) 

A straightforward calculation gives the interaction in the S'-wave channel, 

1/^(k, -k, -k', k') = -2pxA{D/eh) • (8.20) 

Since the interaction is independent of momentum, the k dependence of the superconducting 
order parameter A'^ is small [due to small k dependence of Z{n), see Eq. (|6.3|) ]. The approx- 
imate strength of the dimensionless effective coupling in the singlet channel [see Eq. ( |6.4|) ] 



IS 



,.) _ 2p\A(DleH) 



'-'' - - li^^AWe^ ' ''-''^ 



where we can see immediately that 



-9i'f\<l- (8.22) 



The relation of the transition temperature Tc and the mass enhancement is given by the 
curve g = 2 in Fig. ^ The highest T^ is obtained for mass enhancement about 2.5-3.5. 
Similar calculation for the triplet channel gives 

X sm{kp,^{K,)a)sm{kp^{n')a) + sm{kpy{K,)a)sm{kpy{n')a) , (8.23) 

which results in repulsion. 

B. The case oi U = t = 



The tight binding Hamiltonian (|2.1|) with the choice of parameters corresponding to this 



case has the following form in the k representation: 

35 



H = At'Y, sin(A;,a/2) sm(A;,a/2) (a^,<x%,k,a + h-c.) 

k,(T 



(8.24) 



and can be diagonalized by introducing the 



aa:,k = (pxiX) (C?k,+ + C^k ,-) 

aj/,k = 0?/(k) (c/k,+ - c?k,-) 



(8.25) 



operators, where 0a(k) = sign(sin fcQ,a/2)/v2 and the diagonahzed Hamiltonian is 



H = At'Y,\ MKa/2) sin(A;^a/2)| (4,a,+c?k,a,+ - 4,.,-^k,<x,- 

k,cr 



(8.26) 



With this choice of 05-s, the operators d\^^j^ and rf^ o- + ^^^ associated with the upper band 
which is combined from the upper parts of the two bands labelled by x and y. In the 
following, only the upper band is kept and the index + will be dropped. 

For convenience, the momentum is shifted as k^ ^ k^ — vr/a and ky ^ ky — n/a 
in all of the following formulas. For example, the dispersion relation changes to e = 
At' coskrca/2coskya/2. That shift, however, can be completely incorporated by the am- 
plitudes (psi^n): 



>±x — 9±y 



±i V2 



(8.27) 



in Eqs. {^ and ( pnUQ . 

Introducing the notation Fss = Fq, Fg^s = F2 and F^y 
Fi, the matrix Fgg' takes the form 



-x-y 



-xy 



x-y 



( 



Fi -Fi^ 



F2 Fq — Fi Fi 
Fi — Fi Fq F2 

y — Fi Fi F2 Fq y 
which eigenvectors and eigenvalues are easy to obtain: 



(8.28) 
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=(1) 



.(2) 



1 
1 






,(3) 







1 
1 



= (4) 



(8.29) 



and 



r(i) = ^WT/^ ^ ^(2) ^ ,(2)T/2 ^ r(3) = s(3)^/2 , r^^) = s^'^^/A 



(8.30) 



where the indices d are ordered as {x, —x, y, —y) and the superscript denotes a transposed 
vector. 

The corresponding eigenvalues are 

Ai = U'{F, - F2 + 2Fi) , 

A2 = A3 = f/'(Fo + F2), 

A4 = t/'(Fo -F2- 2Fi) . (8.31) 

The integrals in Fq, Fi and F2 leads to elliptic functions 

p = K'{e^/4t')/irH' 
Fo = (|0.|2> = l/2 
Fi = {iplipy cos{kxa/2) cos{kya/2)) = ep/8t' 



'cosKa) = 1/2 - E'{ep/At')/K'{ep/4:t') 



ing Eq. (f5.21) we get the functions C,'- 


^,{k) = -iV2 


cos{kp^{K,)a/2) + cos{kpy{K,)a/2) 


^2(«^) = -V2sm{kp^{K)a/2) 


^s{f^) = -V2sm{kp^j{K)a/2) 


U{k,) = -iV2 


cos{kp^{n)a/2) — cos{kpy{K)a/2) 



(8.32) 



(8.33) 



and the nonvanishing averages of ^^* are: 
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(|ei(K)P)=4(Fo-F2 + 2Fi)=4Ai 

{\U>^)\^) = 4(^0 -F2- 2Fi) = 4A4 . (8.34) 



The leading terms in Eq. ( p.20| ) with proper symmetry are those with the exponents 

•7(0) 

ijk 

if^ = ilfS = i'/lQ , (8.35) 



2Ai — A2 or Ai + A2 — A4, and the corresponding tjL' are 



or 



rx(0) _ rx(0) _ ry(0) _ ry(0) 



rx(v) _ ,x(u) _ _,y(U) _ _fy('J) _ ^i/tc (q r>c\ 

''124 ~ ''214 ~ ''134 ~ ''314 —''/-'-") {O.OV) 



t' /£)\ (2Ai-A2)p 

r(/^l,AC2,K3) = T7:ei('^l)ei('^2)e2*(^3) - (8.37) 



depending on whether Ep/At > 0.41 or < 0.41 , respectively, so that 

^ei('^i)ei('^2)e2*(^3) (^) 

for Ep/At' < 0.41 and 

2/ / p)\ (Ai+A2-A4)p 

e{m, K2, K,) = - [ei(^l)6('^2) + U^^lMl^2)] a(^3) ( - J (8.38) 

for EpjAt' > 0.41 and for t^^Ki, K2, K3) the C2('«2) should be replaced by ^3(1^2) and multiplied 
by a minus sign for Ep/At' > 0.41. 

The x-s, defined by Eq. (^3|), calculated by using Eqs. (|3.20|) and ( ^.37[ - p73^ ) are k 
dependent with 

p 

Xiii^) = X2(k) = 7^AiA2|6(«^)r 
lb 

for Ep/4t' > 0.41 and 

Xi(k) = W^) = ^A4 [4Ai(|6(A^)r + 16(^^)1') + 4A2|ei(A^)r] 
U^) = j\U>^)\' (8.39) 
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for ep/4t' < 0.41. The average of Xj('^)-s over the Fermi surface is 

^ A2A?/4 , if e^/At' > 0.41 

X = { (8.40) 

A1A2A4/2 , if Ep/if < 0.41 . 

In the case Sp/At' < 0.41, for the singlet superconductivity we get 

V'{^,-^,-k',k') = JlpX,\^,(^)\^\^,(^')\^A{D/\e,\) , (8.41) 

which gives rise to an S'-type superconductivity. The triplet spin channel contribution a 
repulsive interaction of a p-type. 

For the case ep/4t' > 0.41 we get a d-wave repulsive interaction in the singlet channel. (In 



Ref. we got rf-wave attraction, which was due to a sign error. 



IX. DISCUSSION 

The role of the assisted hopping is demonstrated in models where additional to the 
conduction band there are orbitals near the Fermi surface. The occupations of these orbitals 
fluctuate between two values. The state of higher occupation is obtained by adding a heavy 
electron . All the other states are ignored. The energies of these two states include all of 
the intratomic Coulomb interaction, thus these states are fully renormalized in the atomic 
sense. We call the attention to two physical realizations: 

{i) Heavy f -electrons. In this case it must be assumed that one of the renormalized /- 
levels is near to the Fermi energy on the scale of the Fermi energy. If for one of the /-levels 
that condition is satisfled, then the model can be applied by considering the conduction 
electron assisted hybridization of the 4/-electrons with the conduction band. The assisted 
process is induced by the change in the occupation of conduction band in the atomic orbital 
(tight-binding) representation. It is a striking feature, that the very large mass enhancement 
in the conduction band practically eliminates the superconducting state. The moderate mass 
enhancement can be correlated with the superconductivity (see chapter VI.). The possible 
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role of the Coulomb interaction in conduction band has been considered also in Ref. |T3| by 
the slave-boson technique and it makes the superconducting state even more favorable. 

{u)Cu02 plane. The strikingly fiat parts of the electronic band structure are due to 
the non-bonding oxygen orbitals^. E.g. in the YBa2Cu-iOj-5 compound it is either due 
to the apex oxygens taking the tetrahedral positions above or below the Cu site or it is 
generated by the non-bonding vr-orbitals of the oxygens in the Cu02 plane, which are 
oriented perpendicular to the plane in the z-direction. These states can hybridize with the 
orbitals on Cu only if the Cu02 plane is distorted and the Cu and O atoms form separate 
planes. Such a fiat band has been recently observed by experimentEil and also reproduced 
by band structure calculationcHl. The previous case might be related to the influence of 
the distance of the apex oxygen from the Cu02 plane on the superconducting transition 
temperature''. See for more detailed discussion Refs. p!2|-[T^. The calculated induced electron- 



electron interaction is momentum dependent due to the form factors in the electron assisted 



hybridization, but the models treated in Ref. |1^ the transition leads to always s-type of 
superconductors. 

The Hamiltonian for a definite system must be constructed in the atomic orbital picture, 
thus the tight binding approximation is used. For the sake of simplicity only the conduction 
band crossing the Fermi energy is kept for the light particle, even if the other broad bands 
could contribute also to those integrals which have logarithmic character in the most simpler 
approximation. 

The vertex corrections to the conduction electron assisted hopping between the light 
and heavy bands are determined. The dependence on the occupations of the heavy orbitals 
(i.e. interaction like c^K^hh) does not contribute in the logarithmic approximation, thus 
it is not taken into account. The assisted hopping vertices t are strongly renormalized 
by the local Coulomb interaction U between the heavy and light-particles. The general 
formulation is presented in Chapter III., where the number of different couplings are finite 
as the momentum dependence appears only in the form factors which belong to some certain 
class. The solution of these vertex equations for t can be very different depending on the 
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model. In the case of onsite Coulomb interaction U, the renormalization by U can reduce 
the structure in the formfactor by integrating out the dependence on the momenta of the 
scattered electron. That happens in case A in Chapter Vlll. [see Eqs. (8.16)-(8.17)]. On the 
other hand, if the interaction U' is the next neighbor Coulomb interaction, then the form 
factors of the vertex corrections are altering in a certain class, but the level of the structure 
is never reduced, see case B in Chapter VIII. In the one-dimensional model new coupling 
is generated, and as it is discussed in Chapter VII the structure of the vertex equations is 
mapped to those in a one-dimensional interacting electron gas with the two couplings gi, 
g2, gs and (74 known as the g-ologya. The results obtained are generalization of those in Ref. 



TT] and may be relevant in the quasi-one dimensional organic conductors. The functional 
form of the vertex corrections are always power functions of uj/D, but the exponents and 
this way the strengths of the enhancement are very sensitive on the actual structure of the 
formfactors. 

The large mass enhancement described in Chapter IV. and calculated in Chapter VII. 
and VIII. is a quite general consequence of the theory. That can be very large, its value, 
however, is limited by the low energy infrared cutoff due to the dispersion and energy of the 
heavy band. 

The superconductivity is determined by very similar expressions as the mass enhance- 
ment [see e.g. Eqs. ( [4.7|) and ( |6.5|) ]. Usually in other theories the similar expressions are 
related by Ward identities, but here those can not be exploited as the relevant quantities 
appear in different channels (zero sound and Cooper). The single particle weight Z in the 
Green's function are playing crucial role in the strength of the electron-electron interaction 
and the mass enhancement . The related expressions of Chapter VI. and Fig. ^ represents 
quite general relations and they may be relevant in other models as well. 

In the case of one-dimensional models a much richer class of susceptibilities are discussed 
and the phase diagrams in Fig. |T2| contains spin density wave (SDW) and charge density 
wave (CDW) and triplet and singlet superconductivity 

Finally, it is worth to point out, that in the actual calculation the electrons on the 
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different heavy sites are not correlated which is the consequence of our systemic logarithmic 
approximation schema, but that can be lost beyond the approximation applied. As far as 
the number of the excited heavy particle levels at a given time is small, i.e. the dynamics 
occurs on a dilute set of orbitals, the approximation applied is justified. 

For any certain problem with conduction electron assisted hopping between a heavy and 
a light orbital the model can be treated in the general schema presented. The large vertex 
corrections make it promising, that the weaker assisted hoppings can play a determining 
role in some systems, even if their bare amplitude is weaker then the Coulomb interaction^. 
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FIGURES 
FIG. 1. The model for orbitals in one dimensions are shown. The circles represents the light s 

orbitals, and inside them are the heavy d (a) or p (b) orbitals (the / orbitals are not presented as the 

p orbitals have the same odd overlap as the deep / orbitals). The energy levels, the hoppings and 

interactions between them are shown in (c). The clear and shadowed areas indicate the opposite 

signs of the wave functions. 

FIG. 2. The CuO plane of the YBCO compounds is shown as an example of the model in 
two dimensions. The light O and Cu orbitals are found on the sides and on the corners of the 
squares, respectively. The two orbitals of the apex oxygen below the Cu sites can play the role of 
the heavy orbitals. 

FIG. 3. (a) The bare assisted hopping vertex is shown where the double line and light lines 
stand for heavy and light particles, respectively. The wavy line denotes the assisted hopping, (b) 
The bare Coulomb interaction is indicated by dashed lines. 

FIG. 4. (a) The Coulomb corrections are shown in second order. The diagrams contribute by 
logarithms, but they cancel, (b) The corrections appearing in vertex equation are shown by time 
ordered diagrams where the assisted hopping is renormalized by the Coulomb interaction. 

FIG. 5. The contributions to the light particle self-energy are shown by time ordered skeleton 
diagrams. 

FIG. 6. The interaction between the light particles induced by the assisted hopping is shown 
separately for the different channels: (a) the spin parallel and (b) antiparallel channels. The 
diagrams are time ordered. 

FIG. 7. The renormalization of the light electron dispersion curve is shown in the neighborhood 
of the Fermi energy Ep. The renormalization is essential in the range around the Fermi energy 
characterized by the low energy cutoff £h- The large mass enhancement occurs in an energy range 
Zsh. 
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FIG. 8. The selfconsistent equation for the gap A is illustrated. 

FIG. 9. The critical temperature as a function of mass enhancement for different values of q 
[see Eq. Q] 

FIG. 10. The band dispersion of the ID model. 

FIG. 11. The prefactor of (2t cos'^ ip)'^{D /shf'"' of the effective interaction gi, §2 and the x{i^) 
as a function of different filling for the even case: (a) Ai > A3 and (b) Ai < A3. 

FIG. 12. The phase diagram of the ID model in case of (a) even and (b) odd assisted hopping. 
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